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The molar spin susceptibilities X (T) of Na-TCNQ, K-TCNQ and Rb-TCNQ(II) are fit quantita- 
tively to 450 K in terms of half-filled bands of three one-dimensional Hubbard models with extended 
interactions using exact results for finite systems. All three models have bond order wave (BOW) 
and charge density wave (CDW) phases with boundary V = V C (U) for nearest-neighbor interaction 
V and on-site repulsion U. At high T, all three salts have regular stacks of TCNQ~ anion radicals. 
The X (T) fits place Na and K in the CDW phase and Rb(II) in the BOW phase with V « V c . The 
Na and K salts have dimerized stacks at T < Td while Rb(II) has regular stacks at 100K. The x(^) 
analysis extends to dimerized stacks and to dimerization fluctuations in Rb(II). The three models 
yield consistent values of U, V and transfer integrals t for closely related TCNQ - stacks. Model 
parameters based on x(T) are smaller than those from optical data that in turn are considerably 
reduced by electronic polarization from quantum chemical calculation of U, V and t on adjacent 
TCNQ - ions. The x(T) analysis shows that fully relaxed states have reduced model parameters 
compared to optical or vibration spectra of dimerized or regular TCNQ - stacks. 
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I. INTRODUCTION 



The strong 7r-acceptor A = TCNQ (tetracyano- 
quinodimethane) forms an extensive series of ion-radical 
salts [H-Q with closed-shell inorganic ions as well as 
charge-transfer (CT) complexes with 7r-donors such as 
D = TTF (tetrathiafulvalene) . The high conductivity 
and phase transitions of TTF-TCNQ on cooling were 
thoroughly investigated as an important step towards 
the realization of organic superconductivity [4]. TCNQ 
salts crystallize in face-to-face stacks that immediately 
rationalize their quasi-one-dimensional (ID) electronic 
structure. Endres has reviewed the many structural 
motifs of ID stacks |sj. We consider in this paper the 
magnetic properties of "simple" 1:1 alkali-TCNQ salts 
with half-filled stacks of A - radical ions. Complex 
salts with stoichiometry such as 1:2 or 2:3 have less 
than half- filled stacks; they are semiconductors with 
higher conductivity than simple salts. Hubbard and 
related models are the standard approach to TCNQ 
salts or CT complexes [J-Q. Each molecule in a stack 
is a site with a single frontier orbital, the lowest un- 
occupied orbital of A or the highest occupied orbital of D. 

Heisenberg spin chains were initially applied to the 
magnetic properties of TCNQ salts p], especially to 
dimerized stacks whose elementary excitations are 
triplet spin excitons. Subsequently, Hubbard models 
PHI] opened the way to discuss optical and electrical as 
well as magnetic properties. Limited understanding of 
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ID models hampered early treatments. Theoretical and 
numerical advances now make it possible to treat the 
spin susceptibility of ID models almost quantitatively. 
Alkali-TCNQ salts offer the possibility of joint modeling 
of magnetic, optical and vibrational properties. An 
interesting consequence reported below is that different 
model parameters are needed for magnetic and optical 
properties. 

There is considerable literature on K-TCNQ or Na- 
TCNQ, recently in connection with photo- induced phase 
transitions They form g regular stacks with A 

at inversion centers at high T, dimerized stacks with 
two A - per repeat unit at low T. The transitions are 
[9j at Td = 348 K and 395 K, respectively, for Na and 
K-TCNQ. Torrance 110] and others [U sought to model 
Td as a spin-Peierls transition, as discussed in the review 
of Bray et al. 12j who noted that such high Td requires 
unacceptably large exchange constants. The transitions 
have some 3D character since the cations also dimerize 
[HI, [lH. We model the molar spin susceptibility x(T) 
of the Na and K salts at T > Td using regular stacks. 
We also consider x(T) of dimerized stacks for T < Td 
without, however, treating the transition. There are 
two Rb salts: Rb-TCNQ(I) is strongly dimerize d |15l at 
300 K while Rb-TCNQ(II) has regular stacks [lHli3 
with A~ at inversion centers at both 100 and 295 K. 
The recent 100 K structure [13] rules out a dimerization 
transition around 220 K that was inferred from magnetic 
susceptibility [l8| and infrared [lj| data. We reinterpret 
these observations. Regular stacks make Rb-TCNQ(II) 
the best target for modeling x(T). 

Fig. Q] shows the molar spin susceptibilities of Na, K 
and Rb-TCNQ(H). The K and Na data are integrated 
electron spin resonance (esr) of Vegter and Kommandeur 
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[18j, who identified the transitions. Dimerization opens 
a magnetic gap E m > and rationalizes reduced x(T) 
that vanishes at T = 0, whether or not x(^) can be 
modeled. Crystal data QJ, [II] at T > Td indicate 
eclipsed (ring over ring, Fig. [T]) stacks with TCNQ~ 
at inversion centers and interplanar separation R(Na) 
= 3.385 A at 353 K, R(K) = 3.479 A at 413 K. The 
solid line for Rb(II) is esr intensity [llj]. The dotted line 
is static susceptibility [13] corrected for diamagnetism. 
The measurements agree at 300 K and both have a knee 
around 220 K, less prominently in static susceptibility. 
The structure has slipped stacks [17] (ring over external 
bond, Fig. [IJ of TCNQ~ at inversion centers with R 
= 3.174 and 3.241 A at 100 and 295 K, respectively. 
These regular stacks clearly have large E m . They are 
not compatible with finite x(0) at T = and E m = in 
regular Heisenberg [20j or Hubbard [2l| chains. 

The ID extended Hubbard model [Hj] (EHM, Eq. [T]bc- 
low) has nearest-neighbor interaction V in addition to on- 
site repulsion U > 0. Increasing V in a half-filled EHM 
induces a transition to a charge density wave (CDW) 
phase [22]. The CDW boundary is at V c (t = 0) = U/ce M 
in the atomic limit of U » t, where t is electron transfer 
between adjacent sites and qm = 2 for the EHM is the 
Madelung constant of the lattice. The CDW transition 
is closely related to the neutral-ionic transition of CT 
salts from largely neutral DADA stacks to largely ionic 
D+A~D+A~ stacks (23l - [25| . In either case, the ground 
state (gs) undergoes a first-order quantum transition at 
small t/U or a continuous transition when t/U exceeds 
a critical value, or when U < U* . Nakamura (26| rec- 
ognized that the EHM with U < U* has a narrow bond 
order wave (BOW) phase between the CDW phase at 
V C (U) > U/2 and the spin- fluid phase with E m = at 
V S (U) < U/2. The BOW phase has finite E m in a regular 
stack and broken inversion symmetry C; at sites. Sub- 
sequent studies (27l-l29| confirmed the BOW phase of the 
EHM and sought accurate values of V S (U), V C (U) and 
U* . We recently characterized the BOW phase of the 
EHM and related broken Ci symmetry to electronic soli- 
tons [30j. Finite E m in regular stacks is an attractive 
way to rationalize x(T) in Fig. [TJ and we ha ve p roposed 
that Rb-TCNQ(II) is a BOW phase system [Hi Ell. 

In this paper, we model the spin susceptibility of Na, 
K and Rb-TCNQ(II) quantitatively with the EHM and 
related models with more realistic Coulomb interactions. 
We find the Na and K salts to be in the CDW phase with 
V slightly greater than V C (U) and the Rb(II) salt to be 
just on the BOW side of V C {U). Modeling X (T) is both 
challenging and decisive for several reasons. First, the 
full electronic spectrum is required, not just the ground 
state. Second, comparison with experiment is absolute 
since the magnitude of x(T) follows without scaling 
in 7r-radicals with weak spin-orbit coupling. Third, 
all three salts have ID stacks of A - = TCNQ - with 
similar U and other parameters on physical grounds. 
To the best of our knowledge, Hubbard models have 
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FIG. 1: Solid lines: molar magnetic susceptibility of alkali- 
TCNQ salts based on electron spin resonance intensity from 
ref. dotted line, static susceptibility from ref. [13]. 

The Na and K salts dimerize at Td and have regular face- 
to-face stacks for T > Td with separation R between molec- 
ular planes. The Rb salt has regular ring-over-external-bond 
stacks down to 100 K and a knee around 220 K. 



not been applied quantitatively to both magnetic and 
optical/vibronic properties of the same system. 1:1 
alkali-TCNQ salts provide such an opportunity. 

The paper is organized as follows. We present in 
Section II the spin susceptibility of Hubbard-type models 
near the boundary V C (U) between the BOW and CDW 
phases. The magnetic gap E m (V) to the lowest triplet 
state increases rapidly at V ~ V c . In Section III wc 
model the x(^) data in Fig. [T] with similar param- 
eters for TCNQ~ stacks in related but not identical 
crystals. We compute model parameters in Section IV 
for individual TCNQ~ or for adjacent TCNQ~. These 
parameters are reduced substantially in crystals, more so 
for magnetic than for optical or vibrational properties. 
The Discussion briefly addresses the parameters of 
Hiickcl, Hubbard or other semiempirical models. 



II. BOW/CDW BOUNDARY 

We consider a half-filled extended Hubbard model (22[ 
(EHM) in ID and extend it to second-neighbor interac- 
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tions V 2 = jV. The EEHM with 7 > is 

N 
N 

+ Un p (n p - l)/2 
P =i 

N 

+ ^2Vn p (n p+1 +jn p+2 ) (1) 
p=l 

The first term describes electron transfer between 
adjacent sites with retention of spin a. Regular stacks 
in this Section have equal is taken as t = 1. The number 
operator is n p . The last two terms are on-site repulsion 
U > 0, nearest-neighbor interaction V and second- 
neighbor interaction 7V. The spin fluid phase with 
n p = 1 at all sites is the gs for small V while the charge 
density wave (CDW) with two electrons per site on one 
sublattice is the gs for large V. The CDW boundary is 
V c (t = 0) = U/a M (l), where a M (l) = 2(1-7) is the ID 
Madelung constant. As recognized from the beginning 
[H, HH i electrostatic interactions are 3D and ion-radical 
organic salts have qm ~ 1.5. Point charges in ID lead to 
OLyi = 21n2. Physical considerations set olm = 1-5 rather 
than a M (EHM) = 2. 

Finite i in a regular stack leads to a narrow BOW 
phase (2(| between V s < U/au and V c > U/cxm for 
U < U*, withlH u * ~ 7t for thc EHM - Smaller a M 
gives a less cooperative CDW transition and extends 
the BOW phase to higher U* . The point charge model 
(PCM) with long-range Coulomb interactions V n = V/n 
in Eq. ED has 33] C/*(PCM) « lOt. By the same 
analysis, we estimate that the EEHM with 7 = 0.2 
in Eq. [1] has U* ~ 9t. Quantum chemical evaluation 
H of U and V places alkali-TCNQ salts at the CDW 
boundary and imposes the constraint V ks V c in Eq. [1] 

The symmetry properties of ff (7) are the same for 
spin-independent interactions. Total spin S is conserved 
and E m is from the singlet gs to the lowest triplet 
state. The half-filled band has electron-hole symmetry 
J = ±1. We define Ej as the excitation energy to the 
lowest singlet with opposite J from the gs. A regular 
stack has inversion symmetry C, at sites that we label 
as a = ±1 and define E a as excitation to lowest singlet 
with opposite a from the gs. The energy thresholds E m , 
Ej and E a of extended stacks are not known exactly for 
V > in Eq. CD 

We consider V « V C (N) and use valence bond meth- 
ods [HI to solve H(j) exactly for N = An or An + 2 
sites with periodic or antiperiodic boundary conditions, 
respectively. Low-energy excitations are accessible up 
to N = 16, and the full spectrum to N = 10. At 
constant U and t, the condition E m (V) — E a (y) gives 
V S (N) while E a (V) = Ej(V) gives V C (N), We also 



TABLE I: Boundaries V s and V c of the BOW phase of H{y), 
Eq. [T] with TV sites, t = 1, U = 8, 7 = 0.2, and periodic 
(antiperiodic) boundary conditions for TV = 4n(4n + 2) based 
on excitation thresholds E m , E a , Ej. 



N 


Vs(Ecr = E m ) 


Vl{Ea = 0) 


Vc{E a = Ej) 


8 


4.834 


5.105 


5.199 


10 


4.908 


5.139 


5.201 


12 


4.908 


5.139 


5.201 


14 


4.932 


5.149 


5.202 


16 


4.952 


5.157 


5.203 



TABLE II: Magnetic gap E m (V) to the lowest triplet of #(7), 
Eq. [TJ with TV sites, t = 1, U = 8, 7 = 0.2, and periodic 
(antiperiodic) boundary conditions for N — 4n (An + 2) . 



N 


E m {Vi) 


E m (Vc) 


E m (V c + 0.3) 


8 


0.325 


0.496 


1.247 


10 


0.268 


0.451 


1.293 


12 


0.263 


0.413 


1.340 


14 


0.257 


0.384 


1.378 


16 


0.242 


0.364 


1.403 



define Vi(TV) where E a (V) — and the degenerate gs in 
the BOW phase can be explicitly constructed as linear 
combinations of a = ±1 functions 30]. Table U lists 
Vs in units of t for 7 = 0.2 (aM — 1.6) and U — 8t in 
Eq. [1] up to N = 16. The Vs cluster as expected about 
U/cxm — 5.0. Their weak N dependence makes it possi- 
ble to extrapolate to the extended system as discussed 
[36l [37| in connection with a frustrated spin chain. We 
have computed V S (N), Vi(N) and V C (N) of all three 
models (EHM, EEHM, PCM)as functions of U < U* 
and have previously reported [301 EHM values at U — At. 

The magnetic gap E m dominates x(T) as T — > 0. 
It opens [2(| H3| slowly at V s , remains small at Vi 
and grows rapidly on crossing the CDW boundary at 
V c - The size dependence of E m in Table |TT] is for the 
EEHM at U = 8t. Decreasing E m (N) is found in 
spin or Hubbard chains with E m — in the extended 
system. Instead, E m increases with N in all three 
models when V slightly exceeds V C (N). A density 
matrix renormalization group (DMRG) calculation [33j 
for the EHM at U = At shows increasing E m (N, V) 
for N > 30 at V = V C (N). Hence E m of the extended 
system may exceed the N = 8 gaps that we use below for 
V « V c in the BOW phase for V > V c in the CDW phase. 

We compute the full spectrum of if (7) for N = 8 (10) 
sites with periodic (antiperiodic) boundary conditions 
[30I ]. Standard methods give the partition function Q 
and molar spin susceptibility x(T). Fig. [2]shows x( x ) as 
a function of reduced temperature x — ksT/t, where ks 
is the Boltzmann constant. Since #(TCNQ~) « 2.00236, 
the free-electron value, xi x ) i s directly related to Avo- 
gadro's number Na and the Bohr magneton ug. Jiittner 
et al. obtained quantitative x( x ) f° r t ne Hubbard 
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FIG. 2: Molar magnetic susceptibility \ of Hubbard-type 
models, Eq. Q] with finite N. The U = 5t and 6t, V = 
curves for N = 8 and 10 match the extended results (open 
symbols) from ref. [38J] for ksT > 0.3£. The other curves 
have V = V c . The EHM has U = 5t, the EEHM has U = 8t 
and x = 0.2 in Eq. [Tj and the PCM has U = 8t and Coulomb 
interactions V n = V/n. 



model with V = in Eq. Q] and finite x(0); their re- 
sults for U — 5t and 6i are shown by open symbols in 
Fig. [2j The lines are exact N = 8 and 10 results that 
for x > 0.3 coincide with the extended chain within our 
ability to read graphs. The other x( x ) curves in Fig. [2] 
are for N = 8 with periodic boundary conditions and 
V = V c . We find similar x(x) for the EHM with U = 5t, 
the EEHM with U = 8t, 7 = 0.2 and the PCM with 
U = 8t. Small V is conveniently approximated as a 
Hubbard model with an effective U e = U — V. This 
rationalizes reduced x( x ) with increasing V, but not the 
qualitative change of x(0) = due to finite E m in the 
BOW or CDW phase. Large t w 1000X and limited 
thermal stability of ion-radical oganic solids limits x(x) 
to x < 0.5. 



III. MAGNETIC SUSCEPTIBILITY 

In this Section, we model x(T) data in Fig. Q] using 
regular stacks for Rb(II) and for K and Na at T > Tj. 
The first terms of Eq. [T] for dimerized stacks at T < Td 
has transfer integrals 

tp - -(l + (-l) p <5) (2) 

along the stack. We did not change V in dimerized 
stacks. Since all three salts have TCNQ - stacks, similar 
U is expected on physical grounds, and we have sought 
similar U < U* without strictly enforcing the constraint. 



FIG. 3: Three Xm(T) fits shifted by 50 K for clarity of the 
Na-TCNQ data (open symbols) in Fig. [T] for Hubbard-type 
models in Eq. [Tjwith N = 8 and parameters in Table Hill The 
stacks are regular (S = 0) for T > Td, dimerized (S — 0.20) at 
low T and interpolated using Eq. [3] in between. 



It soon became apparent that for T > requires 

V > V c . We studied the PCM with V n = V/n and 
EEHM with 7 = 0.2 in Eq. ED in addition to the EHM 
(7 = 0) in part to search for a fit with V < V c and 
in part to probe the dependence of t and U on the 
model. The following x(T) calculations are all for N = 8 
sites with periodic boundary conditions in Eq. [TJ The 
experimental data in Fig. [T] are now shown as open 
symbols. 

We start with X {T) for Na-TCNQ in Fig. [3] and obtain 
good fits for T > T d for the EHM with U/t = 4, t = 
0.097 eV and V/t = V c + 0.19. The X (T) results for 
the EEHM and PCM are displaced by 50 and 100 K, 
respectively, for clarity. They are equally good for the 
t, U and V parameters listed in Table IIIII All three 
models return t rj 0.10 eV and V slightly larger than V c . 
Good x(T) fits to T = 310 K in the dimerized phase are 
shown in Fig. [3] with 5 = 0.20 in Eq. [5] and the same t, 
U and V. Konno and Saito [l3[ followed the temperature 
dependence of the Na-TCNQ crystal structure and found 
a coexistence region. The regular phase for T > Td = 
345 K appears already at T = 332 K and grows at the 
expense of the dimerized phase that disappears at Td- 
The x(T) fits in Fig. [3] between T x = 310 K and T d are 
linear interpolations according to 

X(T) = |— ^X(T, 5 = 0) + ^-^ X {T, S = 0.2) (3) 
id — ±i la — li 

The coexistence region is 10 K wider in the fit. Ter- 
auchi Q studied the intensity of selected superlattice 
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TABLE III: Parameters for the spin susceptibility of Na, K 
and Rb-TCNQ(II) in Figs. 1513151 



Salt 


Model 


t(eV) 


r r 1 t t\ 

U(eV) 


V(eV) 




T?%4~ A/T 

M/xlM 


U.UaOD 




A Ol A 


Na 


EEHM 


0.0969 


0.630 


0.434 




PCM 


0.0965 


0.627 


0.492 




EHM 


0.0780 


0.370 


0.211 


K 


EEHM 


0.0801 


0.601 


0.406 




PCM 


0.0758 


0.569 


0.440 




EHM 


0.0745 


0.373 


0.199 


Rb 


EEHM 


0.0767 


0.614 


0.399 




PCM 


0.0707 


0.601 


0.440 



reflections for T < T d in both Na and K-TCNQ. The 
intensity is proportional to 5(T) and decreases linearly 
as Td — T near Td- The susceptibility between 320-345 
K can also be modeled as variable S(T). 

Fig. H shows x(T) fits for K-TCNQ, again displaced 
by 50 K for clarity and again in the CDW phase with 
V > V c for T > T d = 398 K. The K-TCNQ parameters 
t, U and V are in Table IIIII The same parameters and 
5 = 0.40 agree with experiment up to 350 K. There 
is no evidence of coexisting phases. The intensity of 
superlattice reflections decreases over an 80 K interval 
and changes discontinuously from 6/2 to at the 
transition @. The solid points are calculated x(Td, 5/2). 
Agreement at x(Td) indicates that x(T) between 350 K 
and Td can be fit with variable 5(T) in these models. 
Smaller t(K) ps 0.08eV^ is consistent with larger R in 
K-TCNQ. 

Figure 5 shows x(T) fits for Rb-TCNQ(II) for the 
three models displaced by 50 K. We took V = V c (8) 
at the upper limit of the BOW phase and set E a = 0. 
The 5 = fit for regular stacks is markedly improved by 
slightly increasing E m beyond E m (8)/t in Table [TTl by 
0.07i for the EHM and by 0.10* for PCM and EEHM. 
Finite-size effects are critical in view of other evidence 
[HIM] for broken Q symmetry in Rb-TCNQ(II), which 
implies V < V c . By contrast, finite-size effects for the 
Na or K salts are absorbed in y > V c . Good 6=0 
fits are obtained down to T « 250 K with the t and U 
parameters in Table IIIII The esr intensity in Fig. [TJ has 
a pronounced knee around Tk n ~ 220 K. The knee is 
less prominent in the static susceptibility. Dimerization 
is ruled out by the 100 K structure, which has regular 
stacks and the 300 K space group |17j . 

An adiabatic (Born-Oppenheimer) approximation for 
the lattice is typically invoked to model the Peierls 39] or 
spin-Peierls [3, S3] instability of ID systems, although 
quantum fluctuations [4l| are important for small 6(0) 
at T = 0. The BOW phase has finite 6(0) for linear 
electron-phonon (e-ph) coupling a to a harmonic lattice 
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FIG. 4: Three Xm(T) fits shifted by 50 K for clarity of the 
K-TCNQ data (open symbols) in Fig. [TJ for Hubbard-type 
models in Eq. [JJwith N — 8 and parameters in Table 1TTT1 The 
stacks are regular (S = 0) for T > Td, dimerized (8 — 0.40) at 
low T and jump to 5/2 at T d . 



[30] |. where a = (dt/du)o is the first term of the Taylor 
expansion of t(R + u). The electronic gs energy per site 
in units of t has a cusp [3fj 

e Q (6)-e (0) = -B(V)\6\+O(6 2 ) (4) 

where B(V) is the order parameter of the BOW phase. 
B(V) w 0.4 is the bond-order difference at V = Vi (17) 
in Table U for all three models for the Us in Table IIIII 
For comparison, a half-filled band of free electrons with 
6 = ±0.1 has comparable B(6) for partial double and 
single bonds. 

The BOW phase has long-range order that cannot 
persist for T > in ID systems. As discussed by Su, 
Schrieffer and Heeger [39[ for free electrons, the extended 
system at low T has regions with reversed 5(0) that are 
separated by topological solitons whose width 2£ goes 
as 1/5(0). Spin solitons are also found numerically in 
the BOW phase [HJ of the EHM or in the magnetic 
properties of organic ion-radical salts [42[ . 

We consider dimerization fluctuations in the BOW 
phase. This regime has equal densities p(T) of spin 
solitons and dimerized segments with successively ±5(0) 
in Eq. [U We approximate each S = 1/2 soliton by 
a regular region of 2£ sites in an otherwise dimerized 
stack. Since E m is not degenerate, E m (V,5(0))/t ini- 
tially increases as B(V)\5(0)\N due to the cusp in Eq. |4j 
as found directly [3(| up to N = 16 at V = V±(N) where 
E a = 0. Such size dependence cannot go on indefinitely. 
It suffices for our purposes to note that 5(0) < 0.10 
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FIG. 5: Three \{T) fits shifted by 50 K for clarity of the Rb- 
TCNQ(II) data (open symbols) in Fig. [T] for Hubbard-type 
models in Eq. [T]with N — 8 and parameters in Table [TTTl The 
stacks are regular (5 = 0). Spin solitons with width 2£ = 60 
in Eq. [5] are used at low T. 



generates large E m with negligible x{T) at low T 
in dimerized regions between solitons. Such regions 
decrease with increasing T and vanish at 2£p(T") = 1 
when the stack is regular everywhere. 

In this approximation, dimerization fluctuations re- 
duce x(T) for T < T. The soliton density p(T) = 
x(T)T/C follows directly from the molar Curie constant 
C = NAg 2 p% °f noninteracting spins. The knee region 
in Fig. [5] up to T is modeled as 



X(T<T',Z) = 2tp(T)x(T) 



(5) 



with 2£ = 60. The fit is adequate for the simple 
treatment of fluctuations. The choice of 2£ = 60 gives 
T" ~ 250 K, somewhat higher than experiment. The 
same soliton width accounts for the T dependence of 
the infrared intensity of a totally symmetric TCNQ - 
vibration [3(| ■ Such IR data is decisive evidence [H, [Zj] 
for broken Q symmetry, whether due to B(8) in well- 
characterized K-TCNQ stacks [H[ at 300 K or to finite 
B(V) in a BOW phase. 

It remains to reconcile dimerization fluctuations at 
low T with the X-ray data for regular stacks at 100 
K and thermal ellipsoids that conservatively limit [46[ 
R + - R_ = 2u < 0.04A. To be detectable, u must 
exceed zero-point motions. The stack at 100 K has small 
p{T) that prevents long range order. Soliton motion 
modulates R as 6(0) = au/t, and the magnitude of a/t 
determines whether <5(0) is consistent with X-ray data. 



We conclude this Section by assessing the parame- 
ters in Table [TTD Three models (EHM, EEHM, PCM) 
with narrow BOW phases have been applied to three 
TCNQ salts (Na, K, Rb(II)). The CT integral t of reg- 
ular stacks depends on overlap, as sketched in Fig. [T] 
and on separation between TCNQ - planes. It is reas- 
suring that the models return identical t to better than 
10 % with t(Na) > t(K) > t(Rb) in an unconstrained 
fit. We sought similar U in TCNQ" stacks. The Us 
in Table UTT1 are identical within 5% for each model. The 
BOW/CDW boundary V c of the EHM with a M = 2 leads 
to U(EHM) < U(EEHM) w U(PCM) « 0.61 eV that 
we prefer on the basis of aM ~ 1.5. The Na and K salts 
are in the CDW phase with V > V c while X (T) of Rb- 
TCNQ(II) is consistent with a BOW phase with V close 
to V c . The x(T) fit of Na leads to 8 = 0.2 up to T « 310 
K and the interpolation in Eq. [3] for coexisting phases up 
to T d . The K-TCNQ fit has larger 8 = 0.4 at low T. The 
knee region of Rb-TCNQ(II) is fit by Eq- [S]with 2£ = 60, 
the soliton width used previously [30( for IR data. 



IV. MODEL PARAMETERS 

The parameters in Table IIIII are for three models with 
BOW and CDW phases. They are internally consistent, 
but considerably smaller than expected from optical 
data. Typical values [H-[l] are t w 0.1 - 0.3 eV and U w 
1 eV in Hubbard models or larger U e = U — V s» 1 eV in 
the EHM. Such parameters rationalize a CT transition 
around 1 eV polarized along the stack and magnetic 
excitations at lower energy At 2 /U e ~ 0.1 eV. We return 
in the Discussion to model parameters. Here we report 
direct evaluation of U, V and t for individual or adjacent 
TCNQ". The results are based on density functional 
theory (B3LYP) with the 6-311**G(p,d) basis in the 
Gaussian 03 package [13] . An eclipsed (TCNQ)^" 2 dimer 
at R = 3.2 or 3.4 A is correctly found to have singlet gs, 
while smaller basis sets [48( yield a triplet gs. Smaller 
basis sets are adequate for model parameters, however, 
as discussed [49[ for t. 

The disproportionation reaction 2A~ —> A 2 ~ + A 
relates U to the gs energies E (A~), E (A) and E (A 2 ~). 
The optimized TCNQ - structure leads to [/(vertical) = 
4.413 eV. Optimization of TCNQ and TCNQ 2- returns 
[/(adiabatic) = 4.192 eV. The relaxation energy of 0.22 
eV for electron transfer is in excellent agreement with 
0.1 eV per TCNQ - deduced [H from Raman and IR 
spectra. The interaction V depends on adjacent TCNQ - 
and can be estimated several ways: (1) electrostatic 
repulsion between the atomic charges qi of the two ions; 

(2) repulsion between obtained in a dimer calculation; 

(3) energy difference 3 E — 2Eq(A~) between the triplet 
gs of the dimer, which precludes the formation of a tt — ir 
bond, and two radical ions. The same values are obtained 
[HI to better than 5%, and Vs in Table ITVl are based on 
the triplet. The listed V(Na) and V(K) are for eclipsed 
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TABLE IV: Calculated model parameters for adjacent 
TCNQ- 



Parameter 


Na-TCNQ 


K-TCNQ 


Rb-TCNQ(II) 


V(eV) 


2.713 


2.671 


2.594 


*(eV) 


0.345 


0.299 


0.182 


h/t 2 (eV) 


0.299 / 0.468 


0.254 / 0.444 




ti/*a (eV) 


0.266 / 0.451 


0.265 / 0.429 




a (eV/A) 


0.59 


0.55 


0.34 



(TCNQ)2- with R = 3.385 and 3.479A, respectively. 
The regular Rb-TCNQ(II) stack has R = 3.241 A and 
a 2.0 A displacement along the long axis shown in Fig. [T] 

The ratio U(ad)/V ~ 1.5 is comparable to au = 1.6 
for 7 = 0.2 in Eq. [TJor to a M {8) = 17/12 for an 8-site 
PCM. The ID stack is close to the CDW boundary of the 
EEHM or PCM. The magnitude of U is strongly reduced 
in the solid state by electronic polarization P m 1 eV 
per charge 50]. Since P is approximately quadratic in 
charge, we have P(A 2 ~) — 2P(A~) w 2 eV. Electronic po- 
larization of adjacent A~ reduces V by a smaller amount. 

The is in Table [TV] are for the 300 K structure of 
Rb and the T > Td structures of Na and K. We find 
t(Na) > t(K) > t(Rb) as expected but for calculated 
ts that exceed the magnetic parameters in Table IIIII by 
a factor of 2.5 for Rb and 3.5 for Na or K. The reason 
for such large reduction is not understood. There are two 
dimerized stacks @ in Na or K-TCNQ at 300 K. Table [TV] 
lists the calculated t\, ti and the larger, smaller separa- 
tion R+,i?_. We obtain 8(Na) = (t 2 -ti) / (h+h) = 0.22 
or 0.26, somewhat larger than 8 = 0.20 at low T in Fig. 
[3) The corresponding 8(K) are 0.27 and 0.25, smaller 
than 8 = 0.40 in Fig. [4] But the K salt has substan- 
tially larger (ti +t2)/2 — 1.1 7i in the dimerized phase 
that leads to an equally good xCO f° r $ — 0-30 when the 
mean value of the transfer integral is used. Overall, the 
calculated and fitted 8 are reasonably consistent. 

The two-point derivative dt/dR — a is an estimate 
for the e-ph coupling constant. The two stacks of Na 
or K-TCNQ at T < Td have almost the same a, whose 
average value is reported in Table [IVJ The 300 and 
100 K structures of Rb-TCNQ(II) return a smaller 
a » O.MeV/A. The structural constraint of 2u < 0.04A 
discussed above leads to <5(0) < au/t = 0.036, an 
estimate that is independent of reduced t in the crystal. 
Dimerization fluctuations of such small amplitude would 
be difficult to detect. 

Direct evaluation of t has been discussed before 
[Hill. Eclipsed TCNQ - gives the largest i(R,0) that 
decreases with increasing separation R. Displacing the 
ion by L along the long axis leads to tilted stacks in Fig. 
[TJ The nodes of the singly occupied orbital of TCNQ - 
generate t(R, L) = at L=1.3A and to secondary 



maxima at other L [49|, [52|. The first maximum at L 
= 2.1 A is close to the Rb-TCNQ(II) or TTF-TCNQ 
structures. A series of substituted perylenes illustrates 
wider variations of t with displacements along both the 
long and short molecular axes [5l"j |. 

We consider next parameters derived from nonmag- 
netic data. Simple TCNQ salts have a broad CT 
absorption Hljct ~ leV polarized along the stack. 
The optical conductivity of K-TCNQ has a shoulder at 
higher energy that has been variously associated with 
dimerization 15311 . with a band edge |54j or with a local 
excited state |4j| of TCNQ". Meneghetti [HJ modeled 
K-TCNQ with special attention to totally symmetric 
mid-IR modes that are coupled in dimerized stacks to 
the CT absorption. Polarized spectra yield the coupling 
constants g n . Meneghetti [55[ used an EHM with A^ = 4 
sites, periodic boundary conditions, and adjustable t\, 
t 2 and Vi, Vi at T < Td- Comparison with experiment 
also entails lifetime or broadening parameters. Nearly 
quantitative fits are shown in Fig. 8 of ref. [55| for 
the optical conductivity at 300 K with coupled mid-IR 
modes and in Fig. 7 for polarized spectra at 27, 300 and 
413 K. 

The EHM parameters of ref. [55] for a regular K- 
TCNQ stack arc t = 0.19 eV, U = 1.20 eV and V = 0.02 
eV. Neglecting V for a moment, we have a Hubbard 
model with x(0) w 1.0 x 10 -4 emu/mole. Including 
V = 0.02 eV in a A^ = 8 calculation leads to x{T) with a 
broad maximum at 1.6 x 10 -4 emu/mole at T max w 800 
K, consistent with the magnitude of regular stacks in 
Fig. [TJ The observed x(T) slope for T > Td is much 
steeper, however, and finite x(0) is not consistent with 
Rb-TCNQ(II). We note that V = 0.02 eV is a finite-size 
effect since A^ = 4 confines an e-h excitation to be close 
together. The CT absorption shifts to lower energy with 
increasing N and optical spectra of longer regular stacks 
return different parameters. The 300 K parameters of 
ref. [55[ are again U = 1.20 eV and alternating i\ =0.10, 
t 2 = 0.37 eV (or 8 = 0.27/0.47 = 0.574) and V x = 0.28, 
V 2 = 0.31 eV. The strong CT absorption at N = 4 
hardly shifts to the red at N = 8 or 12. But 8 = 0.574 
opens a large E m . The calculated x{T) for A^ = 8 with 
these parameters is very small (< 10 -6 emu/mole) up to 
500 K, completely incompatible with the magnetic data 
in Fig. Q] 

Quantitative treatment of e-mv coupling in dimerized 
stacks such as K-TCNQ is based on linear response 
theory and force fields for molecular vibrations (Ha . [57| . 
The coupling constants g n depend on just one electronic 
parameter, the zero-frequency optical conductivity. 
While the CT band is of central importance, its precise 
modeling is not. Dimerized stacks with broken C; 
symmetry are required for coupling to mid-IR modes. 
Coupling to the same mid-IR modes in Rb-TCNQ(II) 
in regular stacks is strong evidence for a BOW phase 
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with broken Q symmetry. The same modes appear 
[HI with slightly higher intensity in powder spectra of 
Rb-TCNQ(I), which is dimerized at 300 K. The T 
dependence of the intensities Im(T) of coupled modes is 
characteristic of a BOW phase, and spin solitons with 
2£ = 60 account ^ for I m (T). 

The optical spectrum in the narrow BOW phase 
is dominated by t due to competition between the 
larger U and V terms. The CT absorption peak is 
around 3t for a regular stack of rigid molecules [3"H 
and shifts to higher energy by [/(vert) - [/(ad) = 0.22 
cV. Dimerization also shifts fiwcT to higher energy. 
Preliminary modeling with all eigenstates of N = 8 or 
10 indicates that the ts in Table IIIII have to be doubled 
for optical spectra and that S » 0.3 produces small 
blue shifts without a shoulder on the high-energy side. 
Larger ts have been assumed all along for optical spectra. 



V. DISCUSSION 

We have modeled the molar spin susceptibility x{T) of 
alkali-TCNQ salts in Fig. [T]that were previously beyond 
quantitative treatment. We have not treated the phase 
transitions of Na or K-TCNQ, but relied on crystal data 
for Td and coexisting phases or diffuse scattering. We 
found consistent parameters in Table UTTl for ID Hubbard 
models with point charges or second-neighbor V that 
reduced the Madelung constant to olm ~ 1-5. The x(T) 
fits in Fig. I3I4I5I have t, U, V parameters in Table IIIII 
that are about half as large as parameters from optical 
data. 

It should perhaps be no surprise that quantitative 
analysis of magnetic and optical data within the same 
model leads to different parameters. Hubbard models 
make the zero-differential-overlap (ZDO) approxima- 
tion of Hiickel theory for conjugated molecules or 
of tight-binding theory in solids. The PCM with 
V n = V/n is a sp ecial case of the Pariser-Parr-Pople 
(PPP) model [11 H3- Salem [go| has summarized 
the merits and limitations of ZDO, which does not 
concern us here. But his discussion of t, the Hiickel 
(3 parameter, bears directly on different magnetic and 
optical parameters. Systematic variations are illustrated 
by many conjugated hydrocarbons with sp 2 hybridized 
C atoms. Hiickel theory provided a convenient approach 
to analyze variations prior to modern digital computers. 
Thermochemical data were successfully fit with a f3 t h 



that is roughly half of f3 op inferred from optical spectra 
foOj . The correlated PPP model with (3 op is defined 
by the geometry of planar conjugated molecules and 
has considerable predictive power [571 l6lj . including 
two-photon spectra and nonlinear optical properties. 
More recently, INDO (intermediate neglect of differential 
overlap) and its spectroscopic version INDO/S have 
different (3 parameters [62j . 

Instead of closely related hydrocarbons, Hubbard 
models are used to study electron-electron correlation 
solids in general. Quantitative application is rare and so 
are homologous series. Moreover, magnetic and optical 
or other properties are typically modeled separately and 
a single half-filled Hubbard band is rarely thought to 
be quantitative. Na, K and Rb-TCNQ(II) are closely 
related quasi- ID systems that nevertheless crystallize in 
different space groups. 

At least qualitatively, differences between magnetic 
and optical parameters may be rationalized in terms 
of relaxed states in thermal equilibrium and electronic 
excitations that are fast compared to atomic or molec- 
ular motions. Equilibrium states that contribute to 
x{T) are fully relaxed with respect to both molecular 
and lattice modes, and relaxed states have reduced 
excitation energies. Hubbard or other approaches to 
electronic excitations start with vertical — excitations. 
Electronic polarization reduces U and V significantly 
in the solid state, but this fast process is fully included 
in model parameters for optical spectra. The Holstein 
model [63| illustrates reduced t due to linear coupling to 
a molecular vibration. Lattice phonons are considered 
in ID for selected modes such as the Peierls mode 
but complete 3D relaxation is prohibitively difficult. 
Yet such relaxation is the most likely explanation for 
small parameters derived from x(T) data. Quantitative 
modeling of the spin susceptibility clearly points to 
different magnetic and optical parameters for Na, K 
and Rb-TCNQ(II). The magnetism also indicates the 
Hubbard-type models for the Na and K salts at T > Td 
are in the CDW phase while the Rb(II) salts is in the 
BOW phase close to the CDW boundary. 
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